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Abstract 

In the usual Fock-and Darwin-formalism with parabolic potential char- 
acterized by the confining energy eg = hujQ « 3.4 meV, but including 
explicitly also the Zeeman coupling between spin and magnetic field, 
we study the combined orbital and spin magnetic properties of quan- 
tum dots in a two-dimensional electron gas with parameters for GaAs, 
for =1 and TV ^ 1 electrons on the dot. 

For N = 1 the magnetization M{T, B) consists of a paramagnetic spin 
contribution and a diamagnetic orbital contribution, which dominate 
in a non-trivial way at low temperature and fields rsp. high tempera- 
ture and fields. 

For ^ I, where orbital and spin effects are intrinsically coupled in a 
subtle way and cannot be separated, we find in a simplified Hartree ap- 
proximation that at iV = m? , i.e. at a half-filled last shell, AI{T, B, N) 
is parallel (antiparallel) to the magnetic field, if temperatures and fields 
are low enough (high enough) , whereas for TV 7^ the magnetization 
oscillates with B and A as a T-dependent periodic function of the 
variable x := 5^^^' "^ith T- independent period Aa; — 1 (where 
m* — 0.067 Too is the small effective mass of GaAs, while toq is the 
electron mass). 
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Correspondingly, by an adiabatic demagnetization process, which should 
only be fast enough with respect to the slow transient time of the mag- 
netic properties of the dot, the temperature of the dot diminishes rsp. 
increases with decreasing magnetic field, and in some cases we obtain 
quite pronounced effects. 

PACS: 75.20.g - Diamagnetism and paramagnetism; 75.30. Sg - Mag- 
netic cooling; 75.90. -|-w - New topics in magnetic properties and ma- 
terials; 73.29.Dx - Electron states in low-dimensional structures 

1 Introduction 

Besides the charge degrees of freedom, the spin of the electrons in 
quantum dots will certainly play an important role in future magneto- 
electronic devices for classical or quantum computing, involving quan- 
tum dots ('artificial atoms'), although the spin degrees of freedom 
are usually neglected, since typically the orbital magnetism dominates 
in quantum dots, as is known, and as we also will see below. How- 
ever, in this paper we look at the magnetic properties of quantum dots 
more in detail, including the 'atypical' spin degrees of freedom, to see 
whether in this way one may be lead to some 'new physics'. More- 
over, it is clear that for our purpose not the most elaborate many-body 
techniques are important, but simple approaches should suffice to draw 
relevant conclusions. With this in mind, we concentrate below on the 
two case = 1 and ^ 1, where A^ is the number of electrons in the 
dot. 

In any case, solids with quantum dots (i.e. planar artificial atoms) 
being placed in an external magnetic field B have to acquire an addi- 
tional magnetic moment. If the dots do not interact, this moment is 
Njy ■ M, where M is defined as the mean magnetic moment of a single 
dot and A^^ the number of dots. That is why the following calculations 
reduce to considering the behaviour of a single dot being in thermody- 
namic equilibrium with the surrounding. In this case we can consider 
the magnetic field B acting on the electrons in the dot as being identical 
to the external field. 

In the following we always assume that the field B is constant in 
space and has the z-direction, whereas the electrons move in the (x, y)- 
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plane. 



2 The case of iV = 1: 

For the beginning, we consider the simple case of a dot with one elec- 
tron. In such a model, and in the usual effective-mass approximation, 
the motion of the electron is described by the Hamiltonian 

n = -^V^ + \m*.\^ + ^ (/; + g*S.) , (1) 

where m* is the effective mass of the electron (= 0.067 ■ tuq for GaAs, 
where mo is the electronic mass), is the Laplacian in two dimen- 
sions, h = h/(27c), with Planck's constant h, ujc = is the cyclotron 
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fequency, uj"^ := a;g -|- where uoq is the parameter characterizing the 
strength of the parabolic potential, which essentially confines the elec- 
tron to the dot. Finally, 1^ = —i{x-^ — y-§^), with integer eigenvalues 
m, is the (reduced) operator of the angular momentum, while the cor- 
responding (reduced) spin-momentum operator 5*2 has the eigenvalues 
±1 ( Here 'reduced' means 'measured in units of W). Furthermore, in 

the following we use the 'effective Bohr magneton' fi*^ := ^ = ; 
g* := — ■ g is the corresponding effective g-iactoi, where for the free 
electron one would have m* = mo and g = 2, whereas for GaAs, we 
have g = —0.44, and (as already mentioned) m* = 0.067 mo. There- 
fore, for GaAs, the quantity \g*\ is -C 1, whereas fi*^ ^ ^b- 
The eigenfunctions of the Hamiltonian (|I]) have the following form, 

V'(f, a) = iPn+,n. {r)Xs. {(^) , (2) 

where Xs^ (c") are normalized eigenfunctions of the spin operator 5*2 with 
eigenvalue Sz = ±|, while the coordinate wave function can be written 

as 

(a)"+ ■ (fe)"- / -|a; + iyp^ 
/oV2vrn+!r2_T V 2/^ 

Here 



^n+,nAr) = , /^^^ . f exp -2 . (3) 



1 

2i 



x + iy _ ^ _^ 
/o ° dx dy 
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with /2 = 



n± =0,1,2, 



±- 



The energy eigenvalues corresponding to these wavefunctions are 



(n: + - + e_ 



(n. 



(4) 



with e± 

The partition function Z{T, B) of the system can be easily calcu- 
lated and is equal to 



Z{T,B) 



E E 

714- =0 n_=0 



ch( 



1/2 

E ^^p( 

=-1/2 

KbT 



— ch(- 



-\ -1 



(5) 



where ks is the Boltzmann constant and T the 'Kelvin temperature'. 
All other characteristic thermodynamic quantities can be found from 
Z(T, i?) by known derivatives. For instance, the mean moment of a 
dot is, II : 

^ = W\{^-^ . (6) 



Thus for AT = 1, the Free Enthalpy G(r, B) = -keT ■ In Z{T, B), 
and the magnetization as well, can simply be separated into 
magnetic' spin contribution, corresponding to the first factor on the 
r.h.s. of d^), and the usual diamagnetic 'orbital' contribution corre- 
sponding to the second factor in (^. Due to the smallness of g* , the 
paramagnetic contribution is very small in GaAs. However at low fields, 
for high enough temperatures the spin contribution dominates in any 
case, since a systematic Taylor expansion shows that the 'paramagnetic 
factor' is Z^pin = 1 + of^ l^^i.ix'^ + •■■ (^•^- correction is oc -B^/T^), 



1 — ^g^4 y4 + ••• (i-e. here the 



whereas the 'orbital factor' is ^orbital 

correction is oc B'^/T^). Note that here we have explicitly used the 

i?-dependence of = \/^^o + which sometimes should not be 'ap- 
proximated away' too early. 
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Of course we are more interested in the low-temperature behaviour: 
In any case, the magnetization M{T, B) can be calculated for = 1 
completely generally from the following formula, with (3 := {ksT)'^, 
and with the characteristic energy cb oc 5, namely := ^y^: 

M(T, B) _ g* (g*eB (3\ ^^^V^" + " 7fe ~ ^^^^^ 
2 ^ V 2 j ch(y'e2 + e|/?)-ch(es/5) 

(7) 

This formula can be evaluated in various limits, i.e. due to the smallness 
of g* for GaAs, one can consider for example the limit ^-y^ ^ 1 while 

at the same time e^/? ^ 1 (i.e. ^ ksT), which is somewhat 

strange, although not unreasonable, if one considers fields in the Tesla 
range and temperatures in the MilliKelvin region. 

In the following, we also consider adiabatic demagnetization or mag- 
netization processes, i.e. where during the change of B and the ensuing 
measuring processes the entropy of the dot is kept constant. This only 
implies that the changes of the i?-field, and the measuring processes 
considered, must be much faster than the thermal relaxation of the dot 
to the surroundings, which is not unreasonable, since with advanced 
techniques magnetic fields can at present be changed significantly in 
two picoseconds, @, |^, whereas the thermal relaxation of the electronic 
state of a quantum dot can be much slower, i.e. by several orders of 
magnitude, 0. 

Now an adiabatic change AB leads to a corresponding change AT, 
which is given by the relation 



(dT\ ^_{g)^^_T_(dM\ ^ 0(dM\ 

Here we have used the well-known relations ^ = ^ (which follows 
from dG = —MdB — SdT) and Cb = (the 'heat capacity of the 
dot' at constant B). So from i.e. from (|^), knowing AB and Cb 
(which must be > and can be calculated from the formula Cb{T, B) = 
UbT ■ 92[TlnZ(T, B)]/dT^), one can directly evaluate AT. 

So at very low temperatures, i.e. if ^^^^ ^ 1 and - of course - 
eo/3 ^ 1, one obtains with the relation th(a;) = 1 — 2 ■ e~^^' + valid 
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asymptotically for x ^ 1: 

dM 1 dM(T.B) , *, . . 

But the heat capacity Cb should remain positive for finite T. So for 
GaAs, in a range of sufficiently low temperatures and sufficiently low 
fields, i.e. for temperatues T below {above) a value Tq{B) given by 



exp 



2to-\g*\eB 



kg -ToiB) 



2 , adiabatic demagnetization (d-B < 0) leads 
to a decrease {increase) of T. If - on the other hand - we do not 
assume \g*\eBP ^ 1, but the opposite limit \g*\eBP <^ 1, then we 
obtain ^ = ■ — ^^]' leading to a similar conclusion, now with 

kBTo{B) = eo • -^^x"' depending on B. 

In Fig. 1, for various values of B, we plot the values of M{T, B) := 
^^'■'^''^^ against the temperature T - ranging from K to ~ 0.008 K 

- and the magnetic induction B - ranging from to 0.08 Tesla; the 
characteristic line Tq{B) separating positive and negative values of M 
is given by the third-lowest contour line from the right, which ends 
for To —> at a value Bk ~ 0.048 Tesla, and for i? ^ at a value 
Tk ^ 0.008 K. 

These are the values for GaAs, calculated with eo = 3.37 meV. {The 
corresponding values for eo = 7.5 meV are: T^ ~ 0.018 K; B^ ~ 0.1 
Tesla, i.e. they scale roughly ~ eo, as expected.) 

In Fig. 2, the adiabatic derivative {^)s from eq. (||) is plotted 
over B rsp. T ranging from to 5 Tesla rsp. from to 6 K in a 3d- 
representation with contour lines. The special contour line separating 
positive and negative values of {^)s does hardly depend on B over 
an extremely wide range of -B-values, and is clearly visible (it is the 
line vertically above the points with T ^ 3 K). In agreement with 
the '3rd principal law of thermodynamics', the adiabatic derivative 
(^)5 always vanishes for T — > 0, for all values of B. But one should 
note that according to Fig. 2 and Fig. 3, {^)s piles up to very high 
values in the region IK < T < 1.8K, for 5- values < 0.016 Tesla: 
Namely, as seen in Fig. 3, in this 'sensitive region' one can easily obtain 
values of the adiabatic derivative between 100 and 500, and even larger 
values for temperatures around 1.5 K, if the external magnetic field 
is around 0.001 Tesla. Note, however, that in our theory we cannot 
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consider naively the limit B ^ 0, since the characteristic magnetic 
lengh, the 'cyclotron radius' lra{B) := ^ ^J^^ = ^ should be much 
smaller than the distance of two dots, or much smaller than any other 
geometrical extension of our 2d GaAs dot system (for B = 1 Tesla, Im 
is 25.7 nm). Keeping this constraint in mind, concerning the change of 
the temperature by an infinitesimal adiabatic demagnetization in the 
above-mentioned 'sensitive region', we have the following result: 

For an isolated qantum dot in 2d-GaAs, with N=l electrons on the 
dot, starting at the point (T ^ 1.5 K, B ^ 0.01 Tesla), for huo 3.37 
meV, we get AT [K] > 100- Ai? [Tesla] . This implies that an unusually 
small adiabatic change of the magnetic field can lead to a significant 
change of the electron temperature in the dot, if one roughly hits the 
above-mentioned region. 

Thus, on the one hand, we have the change of sign of the adiabatic 
derivative from positive values for T < 3 K to negative values for T > 3 
K, a remarkable phenomenon in itself. On the other hand we have the 
fact that the change AT in the 'sensitive region' is unusually large also 
in magnitude, i. e. there it is really important to explicitly consider also 
the spin, and not only the orbital motion. 

Therefore, to diminish the dot temperature (compared with the 
surrounding sohd) e.g. by AT — —0.1 K (at least for a transient time 
tt, which is determined by the small coupling of the dot to the degrees 
of freedom of the surrounding system, and which we assume to be 
much larger than the time tb necessary for significant changes AB of 
the magnetic field), in the above-mentioned region it is only necessary 
to decrease the magnetic field hy B < 10~^ Tesla. 

After having reached the thermodynamic equilibrium of the dot 
with its surroundings, its temperature increases again, but that of the 
solid decreases, until they equalize, i.e. the final temperature has been 
lowered in any case. After that, the magnetic field can be turned off 
isothermally and the process can be periodically iterated. In such a 
way this process can be used for magnetic cooling of the dot system, 
which gives a fiavour of the 'new physics' involved by controlling the 
magnetism of the dot. 

All this will be considerably more effective - and also more inter- 
esting - for a large number of dots and for » 1 electrons per dot : 
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3 The case 1 — a simplified Hartree aproach: 



3.1 Basic approximations 

If the dot contains electrons, the Hamihonian is 

'l '''' 

Here the inclusion of the Coulomb interaction for our multi-electron 
planar parabolic quantum dot leads to considerable complications in 
comparison to free particles, since the Coulomb energy is of the same 
order of magnitude as the kinetic energy for electrons confined in 
dots. This Coulomb interaction is known to consist of the direct term 
(Hartree term) and the exchange term (Fock term): The former in- 
teraction is of long-range type, while the second one is short-ranged 
(cf. |3]-[|1) and oscillatory in it's position-dependence, cf. eq. (13) in 
P|. As a consequence, as shown in a long calculation in |]TU| for which 



we do not have a shorter argument, the ratio of the exchange energy 
divided by the Hartree energy decreases in d=2 dimensions as N~^/^. 
Therefore at sufficiently high N ( > 10^-10^) the exchange should no 
longer play the usual all-important central role - considering also exact 
calculations for C(10) electrons, see e.g. [^j, which show that then the 
above-mentioned ratio is < 1/3. So we neglect the exchange in a kind 
of zeroth-order approximation which still gives interesting analytical 
results for the B- and T-dependence (see below) generalizing directly 
those of the preceding chapter 2. (Including the exchange would pre- 
clude this analysis.) 



Thus, from (10), we arrive at Hartree equations of the form, [11 



(11) 

where q := (r, a); j = 1,2,...,N, and where the lower index pj repre- 
sents a triple of the three quantum numbers n+, n_, and Sz, and where 
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Vjlf) has to be determined self-consistently : 



e J |r — r 

N 

njir) ■■= \^Prir,a)\\ 
i(^i)=i <^ 

e is the dielectric constant of the sohd, e.g. e ^ 12.5 for GaAs. 

(Note that a numerical calculation supports the application of the 
Hartree approach, at least for a qualitative behaviour of multi-electron 
dots, [0].) 

For a solution of the Hartree equations we use as zeroth approxi- 
mation the semi-classical formula for n{r) given in |jl3|, i.e.: 



r2 



n,(f)^n(f)= aS^Vl-^ for r < i? _ (12) 
I else . 



This approximate formula, which admittedly contradicts the boundary 
conditions for the harmonic oscillator functions, describes at least qual- 
itatively the density of electrons inside the dot (comparison with exact 
numerical results for a small number of confined electrons, |]12| , |l^, 
shows that the quantum corrections to n{r) do not modify it essentially 
except near the edge in cases of 'edge reconstruction', as discussed in 
the already mentioned paper 0, see also chapter 4.7 in 0). 

As usual, for parabolic confinement we rely on the close connection 
between the 'effective dot radius i?' used in (|12D and the multi-electron 
wavefunction; therefore one can consider R also as the fitting parameter 
for this wavefunction. 

In zeroth order of perturbation we have with (|T2|) for small r/R : 

V,(r)«-ydVi^« — .^l-^j. (13) 

As a consequence, in the interior of the dot (and not only there, see 
the remark below) we now get for the Hartree differential equation the 
simple 'renormalized form' : 



9 



i.e. through this equation we now have an effective single-particle equa- 
tion, where the 'renormahzed confining frequency' f2, the 'renormahzed 
single-particle energy' Ep, and the 'renormahzed cyclotron length' / are 
defined as 

So all three quantities are now i?-dependent, which should be kept in 
mind. 

Note that for our case, i.e. for ^ 1, one has R ^ I; so almost all 
single-particle wave functions should be exponentially small for r ^ R; 
therefore only a neglegible number of electrons is located near the edge 
of the dot, and the solutions of eq. (0) given by Eqs. @ and (|) can 
also be applied for r > R. 



3.2 Thermodynamics 



At T > the probability to find an electron of the dot in a state with 

37rNe^ ; „ z7 _ . jg defined by the 



an energy 



+ 



Fermi distribution ns(e), i.e. 



1 + exp 



(16) 



where fx = fi{T, N, B) is the chemical potential of the electrons in the 

.1 
■2 



dot, and with s = = ±7 



Us '■= U q — Sz 

p 2 ^2 AeR 



e± := hVt ± hiuJc . (17) 



The chemical potential /i can be determined as usual from the condition 
that N = E iV„ with 



00 00 



Ns{T,B,fi)= J2 ^s(e+n+ + e_n_) 

n_)_=0 n_=0 
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Let us now introduce the Grand Thermodynamic Potentiaiy{T, B, /i) 
E 3^.(T, 5, /i), through 



ys 



E 



—kaT ■ In 



1 + exp 



knT 



with 



N,B 



(19) 



Then the 'mean energy' -E(T, B, N) (i.e. the internal enthalpy^ [0)) of 
the dot), and its Free Enthalpy G(T, B,N) = E — T ■ S, where S is the 
entropy, are determined by the equations 



E 




(20) 



(21) 

Here the final term in eqs. ( pO]) and (|2T|) represents the 'double counting 
correction' of the Coulomb energy, where we have used that 

-.2 , , ^ ^^^^ ^^^/^ ^ 3^ jY2g2 

~ To " 



e' 
Ye 



dV 



d r 



eR 



At T 0, E and G transform into the energy of the ground state of the 



dot: Eq 



lim E - 
lim u: -Rn : = 



limG = 
lim R. 



9tt N^e^ 



eRo 



where 3^o 



hm 3^; 



Now for finite temperatures we define our 'effective dot radius' R = 
R{N, B, T) (or more precise: the 'effective radius of the electron liquid 
on the dot') from the condition that the Free Enthalpy should fulfill 

Again, this condition couples spin and orbital degrees of freedom. 

As in ||T^ , we now use Laplace transforms of the quantities appear- 
ing in (|r^) and (|T^, marked by a 'tilde'; e.g. we write 

c+ioo 



n.Ae 



1 

27ri 



dp ns{p) ■ e~ 



-p-e 



with 
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where c is an arbitrary real number, which must only be 'positive 
enough' to ensure existence of the transformation (see below). Then 
Ns and 3^s can be represented in the following form: 

Ns= I {-^)-Z{e)de, (23) 



ys = -J ri,(e)Z(e)de, (24) 



where the Laplace transform of Z{e) is given by the simple expression 

Zip) = — 7 T . (25) 

Here the constant c in the Laplace transforna (see above) has to be 
chosen in such a way that all singularities of Z{p) are situated to the 
left of the straight line (c — icxD, c + ioo) - this gives a precise meaning to 
the above-mentioned formulation 'positive enough'- : Then the contour 
of integration in (^) can be closed at infinity, and we can use the 
residuum calculus to evaluate the integral. It is easy to show that Z{p) 
has at the same time poles of first order at the points p = p^^ : = 
27rin with n = ±1, ±2, and a pole of third order at p = 0, if the 
quantity X := ^ is ^ Xq, where 



If, on the other hand, the non-generic condition X = Xq is fulfilled, 
then Z{p) has poles of first, second and third order. 

From now on we will be interested only in the 'generic situation', 
i.e. in such fields, for which X ^ Xq. Having found all the residues 
of Z{p) and performing the summation with respect to all poles, we 
find an expression for Z{e). Then, substituting Z{e) into the integrals 



(23) and ( p^ and using the properties of the Fermi functions, we find 
expressions for N and 3^. In the low-temperature limit ksT -C /i<j, 
these quantities take the form 



12 



+ 



+ 



' 2 



E 



PI 



(+) / 



^0 



(-) / f^se+ 



1 



^0 




(27) 



-3 Cq Z 



2^ 



HQ 



4 / huc 



7T 



. /p(+) ^- p(+)(^^s^+\\ 



(-) / /^^e. 



eg 



(2^ 



where := (cUq — 



3tt Ne^ 



X em*R3 1 ' which corresponds to the first equation 
in (p^. Here the periodic functions P^\z) have for even m the form 



2'K^nkTeA 



n=l 



sh 



'2-K^nkTtA 



n -1 



cos(2'7rnz) 



(29) 



whereas for odd m the same resuh apphes, if the functions cos(27rn2;) 
in (p9D are replaced by sin(27mz). 

Finally, we obtain the following expression for the Free Enthalpy : 



(-) / /^^e- 



e- p(+) / 



*~ 2 
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The chemical potential /i is found from (P7|), which can be rewritten as 



1 + 1^ o-fe 




,*\2 



2^ 



TT 

~3 




In the lowest approximation, if the number of electrons in the dot is 
large (A^ ^ 1), then at moderately low magnetic fields {huc ^ 2eo) 
and moderately low temperatures this expression takes the form (in 

zeroth approximation, /i ^ /io): = 1, or /i ^ /iq = eo^/N. 

In the next approximation, within an accuracy of order 1 / ^/N, the 
chemical potential is equal to 



;.-eov^|l + -L/(Ar,i?,r)} 



(31) 



where the quantity 
fiN,B,T) -- 



1 
2 



(9 



*\2 



hujc 

.2^ 



- E 

S=— : 

1 
2 

- E 



pW ^^£l± +p( 



TT 

(-) / ^J^se- 



-P. 



(+) / ^^s^+ 



+ -P2 



(+) / fJ's(^+ 



^0 



is of order-of-magnitude 0{1), and from (|l^) one gets 



4 + 



9 



AeR 
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If we substitute this expression for n{N,B,T) in (pO]), we obtain the 
Free Enthalpy as a function of N , B and T . As follows from (0), the 
first summand in the r.h.s. of (|30D is of order A^2~2o-^ jf ^ jyo- 
can be shown that a = 1/6, see eq. (0) below). The second summand 

All other 



on the r.h.s. of (|30D is of order N"^'"* (with 7 = 
summands have still less order of magnitude, 
first two terms in the r.h.s. of (|30D and use (^ 
following equation for R, |[T0[| : 



3' 



see 



If we retain only the 
then we come to the 




100 a*B 
277iR ^ ' 



(32) 



where a 



B 



is the effective Bohr radius. 



This means that c^o (= ^) is essentially identified with the plasma 
frequency calculated from the electron density in the dot calculated 
with T = and i? = 0; in this approximation the radius R of the dot 
does not depend on B and T at all and is defined only by the number 

of electrons on the dot. To obtain a dependence of R on B and T 
while solving equ. (|3^), it is in principle necessary to take into account 
corrections of higher order; yet at A^ 2> 1 the corrections are very small, 
namely of relative order 0{N~2)^ and we neglect them. 

For a*Q ^ R the solution of (p2[) in the first approximation is 



R^Ro-[l 



100 a*r 



B 



5l7ri?o 



where Rt 







^ 4e m*uj, 







Hence = h- 



N 6, and 



2„2 



20€R 



Ar2 



• 



lOOo^ 
. 27ttR 



)/(l 



• lOOa^ 
. 27ttR 



)) 



oc R 



(33) 



1 

2 oc 



oc oc A^ 3 , as already stated above. 



From Eq. (^) it follows that R transforms into Rq in the classical 
limit {h ^ 0), cf. Eq. (|T2|). Moreover, the dependence of R with respect 
to A^ and Huq as above, corresponds well to numerical results from 
Now the magnetic moment of a dot is defined by the derivative 



M 
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Using expression ( pQ]) for y we find M, taking into account terms of 
order 1/\/N : 



— = Viv- 



(9 



*\2 



* n P*^"*") f w^+ A p(-) w:- A 



1 + 



(34) 

where /io = eo"\/iV has been defined above. 

We see from this expression that M — > for 5^0, since then 
e+ = e_. At the same time the first summand on the r.h.s. of (|34D is 
negative {{g*y < 1) and monotoneously decreasing with increasing B. 
However the second summand oscillates around zero with increasing 



B, and so the possibihty exists that 



2_M 



may be positive at i? — > 0. 



Let us find out conditions when this case can happen: 

For this we expand the r.h.s. of (|3^) into a series in terms of B, and 
in a hnear approximation in B we get for ^ 1 : 



'N 
AN-Po 



*\2 



eo 



+ AN ■ PoiVN) 



eo 



N) 



f^B-B 



eo 



(35) 



The function Po{x) is periodic, with period Ax = 1, namely 

-1 



n=l 



2'K^nkBT 
eo 



sh 



27T^nkBT 
eo 



■ cos(27rna;) 



(36) 



which can also be written as Pn(x) 



+ 00 

E A{x 



n] 



— 1 , where the 



functions A{x — n) are calculated in the Appendix. As a consequence of 
the large- n-behaviour of the Fourier coefficients in front of cos(27rnx) 
in (|SB|), the function A{x) reaches a sharp maximum at x = 0, namely : 



A{x) = 



eo 



AksT 



ch 



f epx 
\2kBT 



A(0) 



eo 



AkBT 



Thus the function Pq{-\/N) in eq. (|35|) takes positive values at A^ 



m 



[m 



0, 1, 2, ...), if at the same time > 1 (see also eq 



below) . 
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Therefore at low fields the quantity is positive, if the following 

two conditions are simultaneously fulfilled, namely (i) the temperature 
has to he low enough: T < Tq := and (ii) the number N of electrons 
in the dot is equal to N = rn^ , with m = 0, 1, 2, ... {which corresponds 
for B = to a 'half filled outer shell' condition, as we shall see). 

In fact, for a planar parabolic dot in the absence of magnetic field 
B and in a one-electron approximation without interaction, the energy 
levels of the electron are defined by quantum numbers n = 0,1,2,..., 
(= n+ +n_ in Eq. (P) and are degenerate with respect to the numbers 



77, 1 — n_ 



0; ±2; ±4; ±n for even n 
±1; ±3; ±5;. ..; ±n for odd 



and, of course, also degenerate with respect to the spin. 

Electronic states of a planar dot with the same quantum number 
n form a 'shell': Then at = the last electron shell (i.e. with the 
highest possible n) of the parabolic dot turns out to be just half-filled. 

(Here it should be noted that in eq. (|35| ) the explicit spin dependence 
- i.e. the term involving (g*)"^ - is neglegible for ^ 1, and the 
phenomenon considered is actually primarily a shell effect.) 

As a consequence, a planar quantum dot with A^ = (^ 1), at 
low temperatures (T < Tq = and at 'sufficiently low fields' (what 
this means quantitatively, is determined soon), i.e. for 5 — > 0, turns out 
to be a 'paramagnetic two-dimensional artificial atom' with magnetic 
moment 

M(T, B, N)=f^l.AN. - l) ■ . (37) 

Thus, when the magnetic field is increased, the magnetic moment of 
a planar dot in a two-dimensional electron gas at first increases oc B 



according to (\3li); then, according to ( \3^ ) it reaches a maximum, then 
diminishes again, vanishes, and becomes negative. 

Let us first find the value Bq, where M vanishes. At A^ = m^, i.e. 



at integer vA^, we have 
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where Pi{x) has been defined above. Then 

|- = V7v.{((9-r-l)^^i:^ + 4ft(V]vM}. (38) 

If T ^ 0, then at < x < 1 : Pi(x) = ^ - x. That is why at T = for 
sufficiently small magnetic fields 

= ViV . ( O* 2 _ ^\ ^ 2 - 4:^/N ■ ^ 

/is ^ eo eo 

The r.h.s. of this expression is positive when 

B<Bo:=^ ^-^ . (39) 

So 'sufficiently small fields' means B ^ Bq; i.e. in the region {B < Bq, 
T < To) the magnetic moment of a planar dot is positive (i.e. the 
dot is paramagnetic) . Outside this region M is < 0, and the dot is 
diamagnetic. 

For GaAs at Huq = 3.37 meV we have the following typical values 
for To and Bq, which should be compared with the results of Figs. 2,3 : 

= 100, To = 3.26 K; Bq = 0.065 T , 

N = 25, To = 4.10 K; Bq = 0.164 T . 



As in the case of a dot with A^ = 1 electron, finally the adiabatic 
temperature derivative (^)g w.r. to changes of the magnetic field is 
calculated through the expression 



dT\ y \^)b,n 

\d-B/^ Cb,n 

where Cb,n is the heat capacity of the dot. From ( ^Tj) and (|^) it is 
easy to show that 

Cb,n = v^^^ + Oi^) . 

-J eo V iV 



If we consider dots with half-filled last electron shell (i.e. A^ = m^), 
then according to ( |5BD 
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(Here one should remember that Pi{x) depends on T, see eq. (^).) 
Hence it follows that 

If we use the relationship ^^^^ = Po{x) and the expression (|36| ) for 
Pq{x), then it is possible to show that at < x < 1 




1 +00 

^i(^) = 2 ^ 



eo ■ (a; - n) 



2kBT 



1 

— x H — 
2 



Hence it follows that the derivative dP{x)/dT vanishes at x=0, |, 1. 
At the same time 

d'Pijx) d'P.jx) eo , 

^ ^ — = — — — — = T < 0, whereas —r—: — > , 

dxdT |x=o dxdT \x=i AksT^ dxdT \x=l 



Thus, at low fields (vA^/i^ ■ i?/eo <C 1), similar as for = 1 (see Fig. 
2), we get positive values of the adiabatic temperature devative, 




eo(/i 



\2 



B) 



27r2 kB ■ (ksTy 



■B>0. (41) 



Altogether this means that the temperature obtained by adiabatic 
demagnetization of the dots, dB < 0, is strongly S-dependent: It first 
diminishes with decreasing B, then reaches a minimum value at i? = 

Bk := eo/(2-\/iV ■ /i^) (since the derivative (^)^ vanishes at ^^b^ = 

i.e. B = Bk), and then begins to raise, due to the inequality (^) < 
at 5 > Bk. 

Concerning the -B-dependence of we mention - however - 

the following point: For A^ = 1, the contour line {^)s = in Fig. 2 
depends on T only, but not on B, for a large region of B values, where 
it is simply given by T ^ 3 K. According to Fig. 2, this seems only to 
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be different for very small 5-values at T > 3 K. In the first-mentioned 
respect there seems to be a qualitative distinction between the cases 
N = 1 and ^ 1, which we do not yet understand at present. 

At last, using expression (PT|), let us find out, how fast the temper- 
ature of a dot diminishes or increases by the adiabatic change of the 
magnetic field : If we let in (^) Ts{B) = T + AT{B), where T is the 
initial temperature, Ts{B) the final temperature, and AT{B) <^ T, 
then it is easy to show that 

AnB).'.(^y.(^y.T. (42) 



47r2 \kBTj V eo 

Thus for quantum dots in a two-dimensional electron gas with GaAs 
parameters at a start temperature of e.g. T = 2 K, and with a 'confining 
energy' eo = hu^ = 3.37 meV, if one wants to change the electron 
temperature on the dot by AT = ±0.1 K for A^ = 100, the B field has 
only to change by ±0.063 T (tesla) - for A^ = 25 by ±0.057 T. If hujo 
= 7.5 meV, then at A^ = 100 the field has to change by ±0.04 T - at 
A^ = 25 by ±0.035 T. 

Finally the following points should be mentioned: as a consequence 
(i) of a pronounced i?-dependence of the total angular-momentum 
quantum numbers L{B) and S{B) of the ground state of the electronic 
ensemble and (ii) of the Coulomb interaction of the electrons, oscilla- 
tions of the physical properties of quantum dots with B have already 
been predicted and observed in a number of papers; e.g. even in an 
early paper of Dingle, Jl^ , before the invention of quantum dots, and 



later-on in papers of Maksym and Chakraborty, 0, and in [|I2|. Fur- 
thermore, numerical results for sma// numbers of electrons, A^ = 0(10), 
show that the maximum electron density may not be at the origin for 
all electron numbers and magnetic fields (see e.g. fl^, [Q, [Q), 



and that the 'electronic edge' of the quantum dot may get a non-trivial 
structure, i.e. the 'edge reconstruction', [|, |^. However, here we stress 
for A^ > 1 that 

(i) in principle the orbital and spin degrees of freedom are intrinsically 
coupled for the individual electrons, although of course the total mo- 
mentum quantum numbers 5* and L remain well-defined for circular 
dots, and that 

(ii) with the quasi-classical electron density (|T^) it follows from eq. (|39| ) 
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that there are not only periodic oscillations of the magnetization of the 
dots with a period oc , but that throughout these oscillations, 
M does not remain negative but alternates periodically in sign. This 
is essentially a 'Hartree shell-effect', and it seems from our analytical 
results that this should be seen parallel to the exchange mechanisms, 
i.e. one is dealing - so to say - with two different sides of one coin. 
This is analogous to the situation in the conventional 3d magnetism, 
where both the 'Hubbard mechanism' (in mean-field approximation es- 
sentially a Hartree effect) and the 'Hund's rule exchange' are important 
for the magnetic properties, although of different relative importance 
for different systems. 

4 Conclusions 

It follows from the above-stated that 

• the magnetization of a planar one-electron dot, = 1, see (P), can 
be separated into two parts: (i) the paramagnetic magnetization 
caused by the own magnetic spin moment of the electron, and (ii) 
the orbital magnetization which is due to the quantized orbital 
motion of the electron of the dot in a magnetic field. At low 
temperatures (T < Tq) and fields {B < Bq) the paramagnetic part 
of the susceptibility exceeds the diamagnetic one, and as a whole 
the dot is paramagnetic whereas at high fields and temperatures 
the dot behaves diamagnetic. An analogous situation is observed 
concerning the adiabatic temperature derivative w.r. to changes of 
the magnetic field: At low fields {B < B^) and temperatures (T < 
Tk), the adiabatic temperature derivative (^)^ is > 0, whereas it 
is < at high fields and temperatures. For more results for GaAs 
parameters. Fig. 2 and Fig. 3 should be consulted, and it should 
be noted that in case of Fig. 3 one obtains quite pronounced effects 
in the 'sensitive region' of T ^ 1.5 K for B < 0.002 Tesla. 

• In a many-electron planar dot (A^ ^ 1), the effects of quantiza- 
tion of orbital motion and the spin effects cannot be separated 
and should be treated simultaneously. We do this within a sim- 
plified Hartree approach leading to a renormalized single-particle 
equation, where the effective radius R of the electron liquid on the 
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dot is used as fitting parameter to minimize the Free Enthalpy of 
the system at finite temperatures. In this way, orbital and spin 
degrees of freedom are now coupled in a rather subtle way. 

With increasing B, for ^ 1, a spin-dependent restructuring 
of energy levels takes place. As a consequence of a shell effect, 
this leads to a periodic change of the magnetic properties of the 
electrons in a dot with varying B, which is a function of the vari- 
able X := ^^|f|^ with period Ax = 1. (Here all parameters have 
their usual meaning, and eo = hujQ is the confinement energy. It 
should be noted that the period does not depend on the tem- 
perature and involves only the characteristic energy scales of the 
system, except from the factor a/TV. The factor itself is of 
course related to the above-mentioned condition of a half-filled 
outer shell, = m^.) Furthermore, the magnetic susceptibility 
of a dot with a half-filled last electron shell changes not only by 
magnitude, but also by sign: For low fields {B < Bq) and tem- 
peratures (T < To) the dot as a whole is paramagnetic, whereas 
with raising B and T the behaviour of the dot transfers from 
paramagnetic to diamagnetic. 

In an analogous manner, the temperature-effect induced by adia- 
batic demagnetization of the dot behaves differently at low tem- 
peratures and fields (T < Tk, B < Bk) and at high fields {B > 
Bk), respectively: In the first-mentioned case, the temperature 
diminishes with an adiabatic decrease of B whereas in the second 
case, it raises. 
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Appendix 

In this appendix we derive the relation between the periodic function 
Po{x), given by the 'Bloch representation' eq. (|36D , and the correspond- 
ing 'Wannier representation' by the functions A{x — n). 
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The periodic function Po{x) (= Pi{x)) is 



Pq{x) = 2 £ ■ cos(2 



+00 



nnx 



n=l 



where 

A{n) 
Using the identity 



2'K^nkBT 



sh 



'2TT^nkBT\ 



i2nnx 



6{x — n) = ^ exp(— i27rria;) 



we find 

Poix) 



+~ +00 

/ ■ e2"'("-") ■ dfc ■ E '^('^ 



-00 

hoc 



+00 



27ri(x— n) 



+00 



1 = A{x-n) - I 



with 



A{x) 



^ksT 



ch 



2A;bT 
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Figure Captions 



Fig. 1: The reduced magnetization M{T,B)/fi*Q is presented as 
a function of the temperature T (in Kelvin units) and the magnetic 
induction B (in Tesla units) for a quantum dot with = 1 electrons 
on it, in a two-dimensional electron gas with the parameters of GaAs, 
and with the confinement potential parameter eo = huo = 3.37 meV. 
Note the change of sign of M from paramagnetic behaviour (M > 
0) to diamagnetic behaviour {M < 0) by crossing the contour line 
where M(T, B) = 0. (The apparent discontinuities of the contour 
lines, representing = -5 • IQ-^, -2.5 • 10-^ ±0, +2.5 • 10-^ 

, as indicated at the margin, are due to inaccuracies of the plotting 
software.) 

Fig. 2: The 'adiabatic temperature derivative' (^)^ is presented 
against the temperature T (in Kelvins) and the magnetic induction 
B (in Teslas) for a quantum dot with = 1 electrons on it, in a 
two-dimensional electron gas with the parameters of GaAs, and with 
the confinement potential parameter Cq = hujQ = 3.37 meV. Note the 
change of sign of the derivative from positive values for low tempera- 
tures (T < 3 K) to negative values for higher temperatures, and note 
the strong increase in the region of 1.5 K for inductions below 1 Tesla. 
(The apparent discontinuities of the contour lines (^)^ — ~1) 

—0.5, ±0, +0.5,..., , as indicated at the margin, are due to inaccuracies 
of the plotting software.) 

Fig. 3: The same as in Fig. 2, but for inductions as low as 0.002 

Tesla and below, where the 'adiabatic temperature derivative' reaches 
extremely high values. (The apparent discontinuities of the contour 

lines (^) ^ (T, B) = 100, 120, 140, , as indicated at the margin, are 

due to inaccuracies of the plotting software.) 



25 



Figure 1 



The reduced magnetization M{T, B) / is presented as a function 
of the temperature T (in Kelvin units) and the magnetic induction 
B (in Tesla units) for a quantum dot with = 1 electrons on it, 
in a two-dimensional electron gas with the parameters of GaAs, and 
with the confinement potential parameter eo = hujQ = 3.37 meV. Note 
the change of sign of M from paramagnetic behaviour (M > 0) to 
diamagnetic behaviour (M < 0) by crossing the contour line, where 
M(T, B) = 0. (The apparent discontinuities of the contour lines, rep- 
resenting ^^^M = -5-10-3, -2.5-10-3, ±0, +2.5-10-3, as indicated 
at the margin, are due to inaccuracies of the plotting software.) 
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Figure 2 : 

The 'adiabatic temperature derivative' (^)^ is presented against the 
temperature T (in Kelvins) and the magnetic induction B (in Teslas) 
for a quantum dot with N = 1 electrons on it, in a two-dimensional 
electron gas with the parameters of GaAs, and with the confinement 
potential parameter eo = huo = 3.37 meV. Note the change of sign of 
the derivative from positive values for low temperatures (T < 3 K) to 
negative values for higher temperatures, and note the strong increase 
in the region of 1.5 K for inductions below 1 Tesla. (The apparent 
discontinuities of the contour lines (^)^ {T, B) = —1, —0.5, ±0, +0.5, 
, as indicated at the margin, are due to inaccuracies of the plotting 
software.) 
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Figure 3 : 

The same as in Fig. 2, but for inductions as low as 0.002 Tesla and 
below, where the 'adiabatic temperature derivative' reaches extremely 
high values. (The apparent discontinuities of the contour lines (^)^ (^7 B) = 
100, 120, 140, , as indicated at the margin, are due to inaccuracies 
of the plotting software.) 
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